Electron-Polarization Coupling in Superconductor-Ferroelectric Superlattices 
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We present a phenomenological model of periodic ferroelectric-superconductor (FE-S) heterostruc- 
tures containing two alternating ferroelectric and superconducting layers. The interaction at the 
FE-S contacts is described as a coupling of the local carrier density of the superconductor with 
the spontaneous ferroelectric polarization near the FE-S interface. We obtain a stable symmetric 
domain-type phase exhibiting a contact-induced polarization and the ferroelectric domain struc- 
ture at temperatures above the bulk ferroelectric transition temperature. With increasing coupling 
energy, we find the appearance of the ferroelectric phase coexisting with the suppressed supercon- 
ductivity in the S-film. The system is analyzed for different thicknesses of the FE- and S-films 
demonstrating the dramatic change of the topology of the phase diagrams with a variation of the 
layers thickness. The results are expected to shed light on processes occurring in high-temperature 
superconducting films grown on perovskite alloy-substrates exhibiting ferroelectric properties at 
lower temperatures. 

PACS numbers: 74.81.-g,74.78.Fk,77.80.-e 
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I. INTRODUCTION 

Epitaxial combination of superconducting films and 
ferroelectric layers is promising for construction of ef- 
fective microwave devices. Especially advantageous is 
the application of small-loss high-temperature cuprate 
superconductors (YBCO)^ and perovskite ferroelectric 
alloys like Ba^Sri^TiOa (BST) and Pb(Zr x Tii_ x )0 3 
(PZT), in heterogeneous miltilayers. As both of 
these compounds have similar perovskite-type crystal 
structure with small lattice mismatch, the well ori- 
ented BST and PZT thin layers can be grown on 
the surfaces of YBCOi^*£*&. In addition, the role of 
strains and off-stoichiometry structural defects at the 
BST(PZT)/YBCO-interfaces can be minimized by the 
use of growing methods such as laser ablation for fabri- 
cation of the miltilayered structures^. 

Since the 1960s, the strong interest in developing elec- 
tric field-effect devices based on ferroelectrics has been 
stimulated by observations of shifts of superconducting 
transition temperature under a switch of the sponta- 
neous polarization. Superconducting field-effect transis- 
tors (SUFETs) cannot be realized with conventional su- 
perconductors since the large coherence lengths and high 
carrier density of the latter result in negligibly small ef- 
fect of the polarization fieldZ^. In contrast, the carrier 
density and coherence lengths of high-T c cuprates (HTS) 
which are two orders of magnitude lower than those in 
conventional superconductors?, have made their use in 
SUFETs very attractive^. 

For example, it was demonstrated in Ref. prill 
that the critical temperature of YBCO films could be 
shifted by about 6 K through the polarization in the 
BaTiC>3 substrate. As was assumed ini2ii!ii£ii&, the in- 
crease of Xjf is connected with the formation of the ac- 
cumulation layers (with charge carrier density enhanced 
due to the polarization directed toward the interface), 



whereas the polarization in the opposite direction in- 
duces a depletion layer with higher resistance and lower 
Tjf . Besides critical temperature, ferroelectric polariza- 
tion affects transport properties of SUFETs leading to 
a change in resistance of 9-25 % for PZT- YBCO and 
PZT-NdBa 2 Cu 3 07-5(NBCO) heterostructuresiZJ£ and 
showing polarization-induced memory effects in resis- 
tance and current, which suggests their possible storage 
applications^. 

In the case when the polarization is parallel to the 
FE-S interface, the use of the zero-P boundary condi- 
tions (P = at the electrode) results in a weak effect 
of the FE-layers scarcely affecting the behavior of the S- 
films. On the other hand, the influence of the FE-layer 
should be essentially strong when the FE-polarization 
vector P is perpendicular to the FE-S interface. Due 
to the specific perovskite-like structure, the properties of 
the high-T c materials are believed to be very close to the 
perovskite-type ferroelectrics. This fact stimulated dis- 
cussions about a coexistence of ferroelectricity and su- 
perconductivity and a possible ferroelectricity of HTS- 
compounds22i2ii2£. As a consequence, the ferroelectric 
polarization has been suggested to be nonzero at the in- 
terface with the HTS-layer with even the possibility of 
penetration inside the superconductor—. 

Recent studies of ferroelectric-HTS layered structures 
were focused on investigations of the role of free P- 
boundary conditions (non-zero P) on dielectric proper- 
ties of FE-layersi2iS. However, little attention has been 
paid to the behavior of superconducting films as well as 
to the effect of the FE-S interface on ferroelectric struc- 
tural transformations. 

In this work, we develop a phenomenological model 
based on the Ginzburg-Landau theory, for the FE-layers 
sandwiched between superconducting films. As com- 
pared to the microscopic approach proposed in Ref. 
where the S-layer was described as two S-planes with a 
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possibility of a charge transfer between the planes and 
the interior of the layer, the proposed model allows one 
to study directly the behavior of the system for different 
thicknesses of the S-layer in the vicinity of the ferroelec- 
tric transition temperature in the FE-film, which is as- 
sumed to be close to the transition temperature of the 
superconductor. We note that for HTS-compounds the 
condition £s(T) ^> £g for the Ginzburg-Landau theory 
is valid over a much wider temperature range than in 
conventional superconductors^. We focus primarily on 
the case when the penetration depth of the polarization- 
induced electric field is larger than the superconducting 
coherence length {Itf/^% ^ 1)> which is also typically 
realized in HTS layers. 

The paper is organized as follows. In section [HI we in- 
troduce the Ginzburg-Landau energy functional for the 
FE-S heterostructure. We derive the Ginzburg-Landau 
equations and elucidate the nature of the boundary con- 
ditions. In section lTTll we analyze in details the results for 
the case when the bulk ferroelectric transition tempera- 
ture lies below the temperature of the superconducting 
transition. We discuss possible phases which can be sta- 
bilized in the system depending on the values of the tem- 
perature and the coupling energy. Corresponding phase 
diagrams are analyzed in section IIVI for different thick- 
nesses of the FE- and S-layers. We discuss here also 
the thicknesses-dependences of the ferroelectric transi- 
tion temperature. Concluding remarks are presented in 
section ??. 



II. THE MODEL 

We consider a system containing periodically alternat- 
ing ferroelectric (FE) layer sandwiched between super- 
conducting (S) films (shown in Fig.^) infinitely extended 
in the x and y directions. We assume here that the po- 
larization P = (0,0, P{z)) is directed perpendicular to 
the FE-S interface. The thicknesses of the FE- and S- 
layers are given by L\ and L2 respectively, and thus the 
total width of the supercell is: L = L\ + L 2 - Generally, 
the whole sandwich structure can be described by the 
following Ginzburg-Landau functional 

F = F FE + F s + F mt . (1) 

The ferroelectric part of the free energy Q is given by 



F, 



FE 



dz 



1, 



a F P 2 (z) + -b F P 4 (z 



1 fdP(z) 
2 CF { — 



(2) 

where a F = a F (T — T C F ) is the inverse of the ferroelectric 
susceptibility and Tjf denotes the bulk Curie tempera- 
ture. For simplicity, we do not include the sixth-order 
term into the P-expansion J2Jl, restricting our analysis 
to second-order phase transformations. Note that since 
the polarization depends on z only, the last term in PJl 

simplifies to VP = dP Jy ■ 




FIG. 1: Scheme of the periodic FE-S multilayer structure. 



For the superconducting S-film without magnetic field 
the free energy can be written as 



F s 



dz 



Li 



\a S mz)\ 2 



1 



-M^{zT + 



2me 



dip(z) 



dz 



(3) 

where as = a s (T — T^) with T c s denoting the bulk super- 
conducting transition temperature, and m, e and |'0(^)| 2 
are the mass, charge and the local density of the Cooper 
pairs respectively. Similarly to (J2J we assume in Q that 
ip(r) depends on z only, i.e., ip = ip(z). Note also that 
without a magnetic field the phase </> = const, and the 
vector potential A = |f ^ = (see the Appendix). The 
Ginzburg-Landau equations in this static case contain 
only the terms with 1^(^)1 & n d below for simplicity we 
use the notations = tj). 

The interface part of the free energy describes the 
interaction between the spontaneous polarization and the 
superconducting charge near the interface. In order to 
derive Fi nt , we consider first the interaction of P with an 
electric field E = (0,0, E(z)) associated with the polar- 
ization distribution^ 



£1 

1 

2 



dz 



1 



PE 



D(z)E(z) 



dzP(z)E(z) + 



47T 
1 

4-7T 



E 2 (z), (4) 



} where D = E + AttP = eE is the electric displacement. 
The field E can be determined from the Maxwell equa- 
tion inside the FE-layer divD = div(E + AnP) = 0. The 
integration in the vicinity of the contacts gives us near 
z = 



E{z) - E(0+) 
and near z = L\ 

E(Li)-E(z) 



-47r(P(z)-P(0+)), (5) 
An(P(L^)-P(z)). (6) 
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In contrast to the normal metal electrodes where P = 
and E — 0, we employ the following two concepts for 
the FE-S interface: (i) P(0~),P(L+) ^ (see the cor- 
responding discussion in Introduction); (ii) an electric 
field can penetrate inside the superconductor (E(z > 
0~), E(z > Lf) 7^ 0) and thus substantially influence the 
superconducting properties^* 1 ^*^. Using these assump- 
tions together with the boundary conditions for DSL 



(D\ s - D\ FE ) ■ n = Aira, 



(7) 



where a denotes the charge density induced by the po- 
larization at the FE-S contact and n is the unit vector 
directed perpendicular to the surface of the ferroelectric 
into the S-layers, we obtain from JSJ and © 



E(Q-) 
E(L+ 
E(z)- 



4ttP(0~) : 
-AtvP(z), 



= -47TO-((T), 

= 47nx(I+), 
+ < z < L± 



(8) 



To derive the expressions for E(z) penetrating inside 
the S-film at a distance of the Thomas-Fermi charge 
screening length Itf, we consider the Poisson's equa- 
tion near the right (Li < z < L\ + Itf) and left 
(-Itf < z < 0) FE-S contacts 2 ! 



,2 d2( P 

l TF^=^ 



(9) 



where — V z ip(z) = E(z), with the boundary conditions 
for the electrostatic potential ip 

- V z <p(Lt) = E(L+), -V^(O-) = E(0"). (10) 

Here l\ F = 6im e i * s l ar g er f° r systems with the lower 
mobile charge carrier concentration no, and ep is the 
local Fermi energy. The solutions obtained from @ , i|l(J|) 
decay exponentially inside the right S-film (z > L\) 

ip(z) = E(L+)l TF ■ E(z) = E(Lf )e"^(ll) 

and the left S-film (z < 0) 

ip(z) = -E(0~)1 T f ■ e 7 ^, E(z) = P(0")e T ^.(12) 

The electric field penetration depth Itf m HTS com- 
pounds is shown to be in the range of about 0.5 nm 
fRef. llOj) . which is significantly larger than Itf ~ 0.02 nm 
in the normal metal electrodes. Note that the total 
charge in the S-layer: Q = J L dzp(z) = includes the 
electron and background contributions: p(z) = Pb(z) + 
p e (z) with the electron charge consisting of the normal 
and superconducting components: p e (z) = p™{z) + p s e (z) 
according to the two-fluid model of superconductivity 2 ^. 
Since we study in this work the superconducting contri- 
bution p s e (z), we assume here that the induced surface 
density a is described by a local increase/decrease of the 
superconducting charge density 



<r(Lt) = Xo 



Li-\-Ijpf 



1^)1^ = 70 !</>(£+) | 2 , 



where 70 ~ XqItf is the coupling coefficient character- 
izing the change of the charge carrier density near the 
contacts due to the spontaneous polarization in the FE- 
layer. 

From the expressions 113(1 . and J3J we obtain 
Pint — F® nt + Ffe-s, where the interaction energy inside 
the FE-layer 



F? nt = 2tt / dzP\z). 



(14) 



The electrostatic interaction of the penetrating field E[z) 
given by (|llfl - H12fl . with the superconducting charge p s e 
near the contacts (Li < z < L x + Itf) an d {—Itf < z < 
0) is determined by 



'FE-S 



dz^{z)\iP{z)\ 2 . (15) 



To calculate the integrals in (|15fl , we use the expansion of 
i/j(z) at a small distance Sz < Itf *C L\ near the contact 
z = L\ 

^z) = ^(Li) + d -^ll{z-L 1 ) + ... (16) 



and z = 



dz 



i>{z) = ^-) + d ^fl Z + 
dz 



(17) 



and the expressions (flip and Ijl2|l for ip. As the result, 
we get 

Fees = -^P(L+)|^(L+)| 2 + Jp(0-)|^(X-)| 2 



dz 



7 



(18) 



3/2 



a(0-)= 7 o|^(^)| 2 , 



(13) 



-<5(|^i+)| 4 + |^(L-)| 4 ) + 2 ( 5 1 (p(Lt)mLt)\ dViL > ' 

-P(O-)WO-)|ffl) + 0( 7 2 ), 

where 7/2 = 47r(l-l/e)^ F , S = 47r(l-l/e)7 ^ F 
and S 1 = 4tt(1 - 2/e)l TF ~ 7 3 / 2 . The first term O 
renormalizes the coefficient ap in and leads to a sym- 
metrical reduction of the deviation of P from its bulk 
value near the boundaries as was shown in Ref. |2^. As 
a first step in the analysis of the role of the interface 
charge carriers, in this work we focus our attention on 
the contribution of the first two interface terms in the 
Ffe-s given by (fTHfl . neglecting the smaller higher order 
terms in the 7-expansion and disregarding the well stud- 
ied term (|14|) . It should be noted that the contribution 
(|18|l can be derived from a microscopic model describing 
an Ising-type FE-layer sandwiched between the S-planes 
with BCS-pairing 24 . In the latter approach, the coupling 
energy — ^\r-r characterizes the change of the one- 
electron energy Ei in the S-plane due to the distortion of 
the nearest atomic group of the FE-layer with the coor- 
dinate R = Rq + (the average distortion becomes 
nonzero in the ferroelectric state). 
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A. The Ginzburg-Landau equations 

Minimization of the total free energy / = J n dil(F F F + 
Fs+FpE-s) where ft denotes the (x, y)-surface, by varia- 
tion of the order parameters P and ip yields the following 
two Ginzburg-Landau equations 



a F P(z) +b F P 3 (z) - c F 



a s ip(z) + b s ip a (z) 



d 2 P(z) 
dz 2 

... h 2 d 2 ip(z) 



me dz 2 



0, < z < Li(19) 
0, Li < z < 1(20) 



The extra surface terms appearing from the variation, 
give the set of the boundary conditions at the right con- 
tcict z — — L-^ 



and at the left contact z = 0~ (z = L~) 



(21) 



c F - 



dP(Q- 



dz 



h 2 #(£-) 

me dz 



+ ji/j(L-)P(Q-) = 0. 



(22) 



Setting 7 = in (|21|l and (|22(l results in ordinary free 
boundary conditions (zero gradients of the order param- 
eters) at the FE-S contacts, which implies in fact the 
continuity of P and ip near the boundaries. However, a 
small finite coupling energy 7^0 gives non-zero gradi- 
ents of P and ip 



dP c ;2 

IWc 

dz 



(23) 



dz 



where ip c and P c denote the corresponding contact values. 

Consider the first expression (|23|l for P. As the right 
hand side is always positive, we conclude that 



dP^ 

dz 



> 0. 



(24) 



For P > this relation gives at the right contact z = L\ 
an increase of P due to the coupling with the surface 
charge, whereas at the left contact z = 0, P decreases. 
Since for P < we have the similar behavior, we conclude 
that independent of the direction of the polarization 



P{z = L 1 ) >P(z = 0). 



(25) 



Let us analyze the second condition (|23[1 for ip which 
in distinction to the latter case depends on the direction 



switching of the polarization (P < 0) results in > 
and in the opposite behavior of the charge density at the 
contacts. 

Below we present the numerical results for the prob- 
lem jn3hl|22J in the case when T C F < T c s , discussing in 
limiting cases the analytical solutions for P and ip. Note 
that the case Tf < T c s is realized in complex insulating 
perovskite alloys such as BST with high Sr content. The 
ferroelectric transition temperature in BST depends on 
the Ba content x, approaching Tf w for x = (SrTiC>3 
is a paraelectric material with no ferroelectric phase tran- 
sition) and Tf = 120° C for x = 1 (BaTi0 3 ). For in- 
stance, for 7.5% Ba, Tf = 60 K 12 which lies already in 
the range of superconducting transition temperatures of 
the HTS compounds. 



III. CASE STUDY: T<f < T<f 

We now present a brief overview of the results ob- 
tained for our model l|19|) - (|22|) for different tempera- 
tures and coupling 7. Fig. [5] shows spatial dependences 
of P and %b normalized by the zero-temperature values 
Pq = y/ci F /b F and ipo = \/ a %/bs for the following 

(2-0 

set of the parameters: t) = rj p % s — 0.2 (77 = 

tb 2 ° 

- LiL - = 0.1 at a low tem- 



a° s = a Q s /a° F ) and 7 = ^ff-^ 
perature r = T/Tf = 0.5, while r c s = T C S /T C F = 1.3. 



Here £g = y 2me ^ii T s is the zero-temperature super- 
conducting coherence length and the correlation length 
for the FE-layer at T = is: £,° F = 



Note 



of the polarization. For P > 0, the gradient < c , m . m u y 

„„a ~- „* + vJ S " la y er ( Fl §- EI b )> lnsct )> V> approaches its bulk value 



that taking into account the possible variation of the 
field penetration depth (section [nj, we can estimate the 
range for the magnitude of the dimensionless coupling 
7: from ~ 10~ 2 (for conventional superconductors) to 
~ 1 (for HTS films). For the films of moderate thick- 
ness (Li/£ F = 1 and £2/^5 = 1), one can observe a 
strong coupling effect demonstrated in Fig. [21 Similarly 
to the theory for bulk ferroelectrics, the two possible di- 
rections of the polarization can be stabilized in the FE- 
layer. As was already discussed in section lll Al (equations 
(|23|l and 11241) ). P increases at the right contact, z = L\, 
and lowers at the left one due to the coupling with the S- 
electrons (Fig. |2fa)). The important effect arising from 
this, is the appearance of two stable solutions for ip in 
the adjacent S-layer (Fig.|2fb)). The first of them (with 
dip/dz < 0, ip{L\) > tp(L)) corresponds to the positive 
branch of P whereas the second one (ip(Li) < ip(L)) 
appears due to the coupling with the negative P. Note 
that the regions with the lower density \ip\ 2 (the depletion 
layer) and the accumulation charge are located near the 
contacts (L — &s < z < L) and (L\ < z < L\ + 5s) 

with 5s 



& = 



h 2 

2meas 



I, whereas inside the thick 



and we obtain an increase of the electronic pair density 
at z = L\ and a lower density at z — L. Furthermore, a ^ b 



as 

bs ■ 
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FIG. 2: Spatial profiles of (a) P and (b) -0 in the ordered 
FE+SC phase at r = 0.5 < t£ = 1 exhibiting the ferroelectric 
ordering in a FE-film and the superconductivity in a S-film of 
moderate thicknesses(Li/^ = 1 and — !)• The insets 

demonstrate the corresponding profiles for thick layers. 



Similar interface effects with an increase of the FE- 
layer thickness are demonstrated in the insets of Fig.[5fa), 
where the regions of the deviations of P from its bulk 

value Pb = \ — t~ are located within a distance 5p ~ 



£f = J\%\ from the interfaces. 

The whole physical behavior becomes highly nontrivial 
for T > Tf (r > 1) in the temperature range t£ < t < 



A. The phase with symmetric domains 

In contrast to the bulk, where P = for T > Tf, the 
coupling with the interface charge leads to rather unusual 
behavior for P and tp as shown in Fig. [21 at r = 1.1. The 
polarization in this case is enhanced, and for Tf < T < 



Tjf and small 7 < 1 is given by 



2c f coshTf^ 



(26) 



where the superconducting pair density at z = L 



m-)\ 2 =Tp 2 b 



1 H — — tanh — - cot — — 

c F b S 2Cf 2£s 



0( 7 4 



(27) 

As follows from l|26(l . in this state a non-zero P ex- 
ists only due to the superconducting electrons, because 
setting 7 = or ip(L) = immediately yields P = 0. 

For the thin FE-films with ip- < 1 or close to the bulk 
transition temperature (where £p - ► 00) , the effect of the 
contacts is particularly pronounced and the polarization 
is linear in z 



P(,) = ^(*-^ML)| 2 . 



(28) 



However, as the FE-layer becomes thicker (|j ^> 1), the 
polarization deviates from zero only near the interfaces, 
whereas inside the layer P — > 0. For example, for z » 



P 



(29) 



and decays exponentially as one moves towards the in- 
terior of the layer. Near the right contact z m L\, the 
behavior is similar 



(30) 



Note that the polarization is antisymmetric with re- 
spect to the center of the FE-layer z — L\]2: P(z) = 
—P(L\ — z), forming two P-domains: with a negative P 
for < z < Li/2, and a positive P for Li/2 < z < L\, 
as shown in Fig. 13a). The center of the layer, where 
P(z = L\/2) = 0, corresponds to the domain wall. 

Consider now the behavior of the S-layer. For thin 
S-films (|a < 1, L 2 + 0) in the vicinity of the T C F 

(Tf < T < Tf), and for a weak coupling 7 < 1, the 
superconducting order parameter is given by 



ip(z) = ip b 



7 



Li 



;-Li-L 2 /2 



2c F b s r] 



tanh-L _£S +0( 7 4). 



26 



sin 



2?s 



, . (31) 

In contrast to the antisymmetric behavior of P in the 
domain-type phase, ip(z) is symmetric about the center 
of the S-layer z = L1+L2/2: ij}{L\+z) = ift(L-z) as can 
be seen in Fig. &b) . We note that for the case of thin 
S- or FE-films < 1 or L 2 /Cs < 1) the effect of 

fluctuations is expected to be essentially important. This 
can be studied by including the additional higher-order 
terms into the free-energy functional ©-© that will be 
a subject of separate work. From 131( 1 it follows imme- 
diately that as 7— > or 77 = — > 00, ip{z) approaches 
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FIG. 3: Spatial profiles of (a) P and (b) ip in the domain 
phase at r = 1.1 > T c exhibiting the weak ferroelecity in 
a thin FE-film and the superconductivity in a thin S-layer 
(Li/^p = 1 and £2/£<j = 1)- The insets show the corre- 
sponding profiles for thick layers. 



its bulk value ipb- Moreover, we obtain the bulk behav- 
ior in |j5T} also when the adjacent FE-film becomes thin, 
so that L\ <C £f- Similarly to the FE-layer, in a thick 
S-layer the influence of the interfaces is remarkable only 
near the contacts (Fig.|3fb), inset). 

To gain a better understanding of the nature of the 
domain phase, we determine its free energy which can be 
obtained by the substitution of P and ip from l|26[) and 
(ETJ into ©, © and JTSJ 



7 

4c F 



j tanh 



2£ F 



T bs /4 



0( 7 4 ). (32) 



Note that since for the free (non-coupled) FE-layer 
F FE = for T > Tjf , the total free energy for the non- 
coupled FE+S multilayers for T C F < T < T c s is given by 



r FE+S 



T b S ,4 



(33) 



Therefore, the domain state always has an energy gain 
due to a negative difference between the energies in these 
two phases 



^Pdrnn F 



dmn ' 



~ r FE+S 



£f?/> 4 tanh — — • 
4cf 2£ F 



-0( 7 4 ). 

The latter result supports our conclusion that a non-zero 
coupling with ip ^ is the main reason for the stabiliza- 
tion of the domain state. We note also that as 7 — > 0, we 
obtain that the energy gain AF^mn — > 0, the polarization 
P — > 0, and hence the domain phase transforms in fact 
into an ordinary paraelectric state: 



P 



F,7 



T 



as 7 



0. 



Thus in the following analysis we define Tf n as the tem- 
perature of the transition between the ordered FE+SC 
phase with the bulk-type F / 0, and the domain-type 
phase, which occurs on heating. 



B. Strong coupling behavior 

As results from the previous section, generally a weak 
coupling 7 -C 1 at the contacts produces symmetric do- 
main phase with a non-zero polarization for Tf < T < 
Tjf . We demonstrate how this behavior changes as the 
coupling becomes stronger, i.e., with the increasing 7. 
In Fig. 21 we plot the order parameters vs 7 for layers 
of thicknesses L\/£ F — 2 and ^2/^5 = 2. At lower r 
(t = 0.5 and r = 0.55), the behavior of ip and P rapidly 
changes at 7 rj 0.33, corresponding to the transition to 
the normal (non-superconducting) state (ip = 0). In this 
state, P = Pb = const does not depend on z because set- 
ting ip = in (|22|) and (|21fl gives the free boundary con- 
ditions. The suppression of the superconductivity due to 
the strong contact with polarization resembles the effect 
of an external field, impurities and other pair-breaking 
factors on superconductors^. Like these factors, the po- 
larization at the boundaries with the S-layer supplies ad- 
ditional energy l|18|) for the breaking of electron pairs in 
the strong-coupling regime. 

From the point of view of the Landau- Ginzburg func- 
tional ©, J3J), CBJi, 

an increase of 7 induces larger val- 
ues of the fourth-order fluctuation terms as well as the 
third terms in J2J and © describing the spatial inho- 
mogeneities, and hence drives a transition to the homo- 
geneous state with ip = which has the lower energy. 
To explain qualitatively the effect of 7, we consider for 
simplicity only the left FE-S contact with the energy 
TjP(0~)ip 2 (L~) assuming at the right contact the free 
boundary conditions. At low T the expression for the 
free energy in the ordered FE+S-phase is given by 



Ffe 



+s 



■iPtbsL 2 ] + jP b iPl 



(35) 



The critical value for 7 can be found by equating l|35|) 
and the free energy for the state with the suppressed 



7 



superconductivity 



suppr 



1 

T 



which yields 



7c 



(36) 



(37) 



From Ij37fl it follows, that as the S-layer becomes thicker, 
a stronger critical coupling is required to destroy the su- 
perconducting state. In contrast, a larger polarization 
P in Eq. (|37|l reduces which supports our interpreta- 
tion of the analogy with other external field pair breaking 
effects. 

Consider now the 7-dependence at higher tempera- 
tures shown in the insets of Fig. Ufa) and Fig. Efb) 
slightly below (r = 0.9) and above (r = 1.15, where 
the system for small 7 is in the domain phase) the bulk 
value Tjf = 1. Like at low temperatures, the supercon- 
ductivity of the domain phase (t = 1.15) is destroyed 
when approaching the critical value 7^ ks 0.55. We note 
that since i/j = for 7 > 7®, the polarization, which is 
enhanced here only due to the coupling with ip, vanishes 
and we obtain the combined paraelectric+suppressed su- 
perconducting (P = 0, ip = 0) phase for r > . We 
emphasize that in contrast to the bulk ferroelectrics ex- 
hibiting P = in this temperature range, in our coupled 
FE-layer the zero polarization is only possible due to the 
strong coupling effect depressing the source for the en- 
hancement of P. 

The behavior for large 7 can be even more compli- 
cated as demonstrated in the inset of Fig. 0] at the 
temperature r = 0.9 which is below but close to 
(Fig. EI a ))- Despite r < ', the 7-dependence of P 
and tp clearly shows the stabilization of the domain state 
with P(0) = -P(L\) and ip(Li) = ip(L) in the range 
(0.34 < 7 < 0.46), and only with the further 7 increase 
(7 > 0.46), the transition to the (P = P&, ijj = 0)-state. 
To analyze more comprehensively the behavior for differ- 
ent 7, we discuss in the next section the phase diagrams 
(7, t) plotted for the different FE- and S-layer thick- 
nesses. 



IV. PHASE DIAGRAMS (7, r) 

The phase diagrams (7, r) plotted in Fig. [5] and Fig.0 
show the regions of stability for different phases discussed 
in section IIIII As has been revealed in the foregoing 
analysis, the following phases can exist in our periodi- 
cal multi-layer system: (i) the ordered bulk-type FE+SC 
phase with the nonzero spontaneous polarization in the 
FE-layer and the superconductivity in the S-layer; (ii) 
the superconductivity in the S-layer coexisting with the 
ferroelectric domain state exhibiting the symmetric (neg- 
ative and positive) polarization domains due to the FE-S 
contacts (so called domain-type phase); (iii) the phase 
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FIG. 4: (a) Polarization in the FE-film of thickness Li/£° = 2 
and (b) tp at z = L in the S-film of thickness L2 / £,% = 2 the 
magnitude of the coupling 7, exhibiting a transition from the 
FE+SC phase to the FE-state with suppressed superconduc- 
tivity. Here fj — 0.5. The insets show the P and ip behavior at 
the FE-S contacts with increasing 7 for r = 0.9 (solid curves 
and dashes) and r = 1.15 (open squares) with the stabiliza- 
tion of the domain state for 7 < 7*. 



with the ferroelectricity and the suppressed superconduc- 
tivity (FE+NM); (iv) the high-temperature paraelectric 
(in the FE-layer) +normal state (in the S-layer). 

In Fig. [3] we plot a set of the diagrams for the S-film 
of moderate thickness £2/^5 = 1 coupled to the FE- 
layers of different thickness: (a) a thin FE-film, L\/^% = 
0.3,1.0, (b) a FE-layer of moderate thickness Li/£p = 
1.5, and (c) a thick FE-layer with = 2.5,5.0. We 

associate in the following analysis the variation of the FE- 
thickness as related to that of the S-layer, to the ratio 



c 



Li 
L 2 



L2/e s ' 



(38) 



where the ratio between the superconductor and the fer- 
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roelectric coherence lengths 77 ranges from ~ 0.05 (for 
HTS compounds) to ~ 10 (conventional superconduc- 
tors), taking the typical value for £p w 10 nm fRef. Efflh. 
As C decreases from 5/r] (Fig. Etc)) to O.3/77 (Fig. Eta)), 
the area of the domain-type phase drastically expands to- 
wards larger 7 and lower T and finally approaches even 
the low-temperature range as represented by arrows in 
Fig. Et a ) ■ The expansion of the domain phase is accom- 
panied by the simultaneous decrease of the region of the 
FE+NM(suppressed)-phase. We note that in the system 
with a thicker FE-layer (L±/£p = 1, Fig.Et a )), the tran- 
sition to the suppressed state occurs at a significantly 
lower 7^ (7® w 0.9) as compared to a thin FE-film (case 
= 0.3 in Fig. Eta)) with 7^ w 1.8 (r = 1.1), which 
clearly demonstrates the crucial role of the ferroelectric 
layer in the suppression of the superconductivity. These 
modifications of the phase diagram can be explained by 
the competition between the first tendency to destroy 
the superconductivity prevailing for the coupling with 
thicker FE-layers, and the second tendency to retain the 
superconducting state and stabilize P-domains dominat- 
ing in systems with thicker S-layers. As is demonstrated 
in Fig. El the decrease of £ results in a dominant role of 
the S-layer and, as a consequence, in the enhancement 
of the region with tp ^ at lower T < Tf (which in 
turn induces P-domains in the FE-layer). Since in the 
HTS-films r) is about three orders of magnitude smaller 
than in conventional superconductors, one might expect 
from i|38|) that even a slight variation of the S- or FE-film 
thickness leads to a dramatic modification of the phase 
diagram and, consequently, of the stable state of this type 
of compounds. 

Fig. El shows the behavior of the system with a fixed 
FE-thickness Lx/£p = 1 and with a S-thickness i^/Cs 
varying from 15 (Fig. Eta)) to 0.2 (Fig. Etb)). As in the 
previous case, we find here the same tendency of expan- 
sion of the domain phase with decreasing £. We also bring 
to attention that the destruction of the superconducting 
state demonstrated here has been observed in the mea- 
surements of the resistivity in the film YBCO/BaTiC>3 
composites^. 



Comparison with the microscopic model of 
Ref.El 



As has already been mentioned in section |E] when the 
ferroelectric critical temperature Tf is close to T c s , the 
phenomenological approach proposed in this work, has 
common key as pec ts with the microscopic description 
proposed in Ref. .2.4 However, while the latter model de- 
scribes the superconducting layer as two boundary 2-D 
planes exhibiting BCS-superconductivity, the present ap- 
proach allows us to study the behavior inside the thick S- 
layers. Despite the fact that the microscopic BCS-model 
of Ref. does not consider directly the internal pro- 
cesses in the superconducting layer and treats the inner 
part of the layer as a source for the electron transfer to- 
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FIG. 5: Phase diagrams (7, r) for a compound containing 
S-films of thickness L^/^s = 1.0 coupled to the FE-layers 
of thickness: (a) L x /^% = 0.3 and Li/£& = 1.0 (thin FE- 
films); (b) Li/$p = 1.5 (moderately thick FE-film) and (c) 
= 2 -5. 5.0 (thick FE-layer). Here 77 = 0.5. 
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B. Dependence of T C F ' 7 on the FE- and S-layer 
thicknesses 

The behavior of the different phases discussed above 
with the variation of the FE- and S-thicknesses results in 
different types (from monotonous to complex nonmono- 
tonic) of the corresponding T F ' 7 dependences. 

Since the system properties drastically change for the 
large coupling 7 > 7-, the L\- and Independence of T C F ' 7 
is also qualitatively different in the weak (7 < -%) and in 
the strong coupling regime (7 > %). To illustrate this, 
in Fig. ?? and Fig. ?? we plot several curves for various 
values of 7. Specifically, Fig. ?? (inset) and Fig. ??(a) 
show the case 7 < %, whereas 7 > 7 C is represented by 
Fig. ??(7 > 0.9) and Fig. ??(b). 

Depending on the value of 7, we can obtain the follow- 
ing characteristic types of T F ' 7 (Li) and T F7 (L 2 ) behav- 
ior: 

(1) for 7 < 7c, rf' 7 = T F7 /T F decays monotonically 
with increasing L\ to the bulk value r F = 1 (Fig. ??, 
7 = O.f in the in set) . We note that in contrast to 
the model in Ref. :2J where T F ' 7 lies slightly below 
rpF,buik > th e coupling with a S-layer of a finite thick- 
ness produces here the enhancement of the bulk-type 
FE+SC-phase. This corresponds to the increase of T F7 
(Tf < Tf- 7 < T c s ) for 7 ^ (7 < 7 C ), observed also in 
the phase diagrams in Figs. and As a result, we also 
obtain that: 

(2) for 7 < 7c, the increase of Li produces a monotonic 
increase of t F7 , with finally approaching a saturation 
value which is the higher the larger 7 (Fig. ??(a)). 

(3) For 7 > 7c, the complex behavior of r F ' 7 vs 



FIG. 6: Phase diagrams (7, r) for a compound containing FE- 
films of thickness Li/^p = 1.0 coupled to S-layers of thickness: 
(a) L 2 /£% = 5.0 and L 2 /£°s = 15.0 (thick S-layer) and (b) 
L 2 /i% = 0.2 (thin S-film). Here rj = 0.5. 



wards (away from) the boundary S-planes, the following 
key results appear in both approaches: 

(i) the enhancement of the ferroelectric polarization with 
a formation of positive and negative P-domains above the 
ferroelectric transition temperature as well as for large 7 
and low temperatures; 

(ii) the suppression of the superconductivity for the 
strong FE-S coupling due to the ferroelectric field-type 
effect. 

Thus, the basic structure of the phase diagrams (7, T) 
together with the behavior of the system in different ob- 
served phases discussed in Ref. |2J, also retains its fun- 
damental features in the present phenomenological ap- 
proach. 
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FIG. 7: Ferroelectric transition temperature vs the FE-layer 
thickness (for a S-film with I/2/Cs — 1-0) for f] = 0.5 in the 
strong coupling regime. The inset shows r F ' 7 in the weak- 
coupling regime for 7 = 0.1. 

L\ shown in Fig. ??(7 > 0.9), is governed by a de- 
crease of the domain phase-region as Li increases (see 
the plots in Fig. [SJ. First, for thin FE-films, the sys- 
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similar, but reversed "crossover" between the two above- 
mentioned different transition types. Specifically, while 
for thin S-films, the increase of T results in the transition 
from the FE+NM(suppressed) to the paraelectric state 
(at Tjf' 7 = 1), the increase of the S-thickncss induces 
the enhancement of the domain-type phase, leading to 
a lower r^' 7 (which now corresponds to the transition 
between the FE+NM and the domain phases). See the 
evolution of the phase diagram topology as L 2 /^ Q S varies, 
shown in Fig.EJ Furthermore, we note that for the larger 
L 2 /£s > 0.7 the transition from the FE+NM- to the sym- 
metric domain-state exhibits: 

(i) monotonic decay of t^' 1 and vanishing at L 2 /t^ s = 0.7 
(case 7 = 0.67 in Fig. ??(b)); 

(ii) reentrant behavior characterized by the vanishing of 
rf' 7 in a certain interval £ 2 /£<j (1 < L 2 /£,% < 8.2 for 
7 = 0.75 in Fig. ??(b)) and finite r^' 7 otherwise; 

(iii) for the larger values of 7 (7 = 0.9 in Fig. ??(b)), 
77^' 7 first decays to a finite value exhibiting a minimum 
at a particular L 2 (L 2 /^g « 3), and then increases again 
with a further increase of the S-thickness, approaching 
a saturation value for thick S-layers. This behavior has 
a simple physical explanation: once the domain state is 
expanded down to T — (yielding r^ 7 = at 7 = 0.75 
and L 2 /£g = 1.0) the further increase of L 2 for strong 
coupling leads to the larger spatial inhomogeneous terms 
in the FE-energy @, which results in a slight enhance- 
ment of the suppressed state towards higher T and as a 
result, to an increase of t^ 1 (see the two phase diagrams 
for L 2 ji% = 5 and L 2 /£° s = 15 shown in Fig. Eta)). 



CONCLUSIONS 



FIG. 8: Ferroelectric transition temperature vs the S-layer 
thickness (for a FE-film with = 1.0) for fj — 0.5 and dif- 

ferent values of coupling strength 7: (a) weak coupling regime 
and (b) strong coupling regime. 



tern is in the domain phase already for T ~ 0. As 
Li/Cp increases, the area of the domain state shrinks 
towards higher T and finally disappears, which results 
in an increase of r^' 7 up to rf = 1 related to the 
transition from the FE+NM (suppressed) to the paraelec- 
tric state (see the corresponding diagram in Fig. E^c)). 
Note that for a small thickness of the FE-film (Lx/£p ~ 
0.5 — 1 in Fig. ??), the behavior of rf n reflects a 
" crossover" between these two different ferroelectric tran- 
sition types, namely, between the transition from the 
FE+NM(suppressed) to the domain state for the thin 
FE-films, and the transition from the FE+NM to the 
paraelectric phase for the thick FE-layers. This effect is 
also seen in the large 7-region of the diagrams plotted in 
Fig.0 

(4) Fig. ??(b) shows a nonmonotonic dependence of T ( f' 7 
on L 2 for 7 > As in the Li-dependences, for a given 
value of 7 at a moderate thickness £i/£s ~ 1, we find 



In this paper we have developed a phenomenological 
model describing periodic ferroelectric-superconductor 
multilayers. This model generalizes in some aspects the 
microscopic approach proposed in Ref. |3 where the su- 
perconducting layers are represented by two planes with a 
possibility of charge transfer between the planes and the 
interior of the layer. We analyze here the case when the 
ferroelectric transition temperature lies below the tem- 
perature of the transition to the superconducting state 
(T C F < Tjf) which is realized in some complex insulating 
perovskite alloys such as Ba^Si'i-^ Ti0 3 (BST) coupled 
to high-temperature superconductors. Similarly to the 
microscopic model, we observe the stabilization of the 
symmetric domain-type phase exhibiting a weak polar- 
ization due to the FE-S contacts above the bulk ferro- 
electric transition temperature. On the other hand, we 
find the polarization-induced suppression of the super- 
conductivity in the S-layer. We show that the ultimate 
reason for the destruction of the superconducting state 
for strong FE-S coupling, is the increase of the contribu- 
tion of the spatial inhomogeneities to the energy of the 
system. 

We study the behavior for different thicknesses of the 
FE and S-layer demonstrating the strong change of the 



11 



topology of the phase diagrams (7, r) with a variation of 
L t /e F and L 2 /e s - 

In general, we find several types of T ( f ,7 (Li) and 
Tf n (L2) behavior. For the weak FE-S coupling, 7 < j c , 
we obtain a monotonic variation of X^' 7 , with approach- 
ing saturation value in the system with the thick FE or S- 
layer. In contrast to this, for the strong coupling, 7 > 
the dependence of T C F ' 7 vs layer thicknesses appears to 
be highly nonmonotonic due to the competition between 
the first tendency to expand the domain state (prevail- 
ing for the thick S-layers) and the second tendency to 
destroy the superconductivity dominating for the thicker 
FE-layers. 
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APPENDIX: DETERMINATION OF PHASE <f> IN 
THE CASE WITHOUT MAGNETIC FIELD 

Let us consider in more details the superconducting 
layer without a magnetic field. The general form of the 
Ginzburg-Landau functional in the S-layer can be written 
as 



+ 



2me 



dz 

d 

dz 



las^z)? + \b s \i>{z)\ A 



2ie 



(39) 



he 



A{z) i>{z 



Variation of (??) with respect to A, if> and ip* gives the 
following Ginzburg-Landau equations for the determina- 
tion of A, and <fi: 



(40) 



^A(z)-<f>'(z)=0, 
he 

h 2 

(as + bs\ip\ 2 )\ip\ COS0 

me 

x (tp" cos 4> - ip'<f>' sin <j> - (4>') 2 ip cos 4>) = 0, 

h 2 

{as + bs\ip\ sin0 

me 

x {ij}" sin (j> + ip'<f>' cos <f> - (<j)') 2 ip sin <p) = 0. 



Note that in the obtained in this case boundary con- 
ditions the terms containing <fi and </>' appear with the 
opposite signs and thus cancel, resulting in the standard 
form given by Ij21|l and l|22(l in the paper. 

From (??) we immediately obtain that the single pos- 
sible solution for <j> is: 

4> = 4>q = const, L± < z < L. (41) 
We assume here for simplicity the absence of the tunnel- 
ing current through the insulating FE-layer, which cor- 
responds to the zero phase difference on the contacts be- 
tween the S-films: 



= <f>(L+) - 0(0-) = 0, 



(42) 



and leads to the constant phase in the whole system: 
<j> = <f>Q. As the result, we get from (??) 



A(z) = ^4,'(z) = 0, 
le 



(43) 



where ijj(z) = \ip(z)\exp(i(j)(z)). 



and the next two equations transform into the equation 
(PU|l for the absolute value of the superconducting order 
parameter |-0(z)|. 
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